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The discovery or exclusion of the fundamental standard scalar is a hot topic, given the data of LEP, 
the Tevatron and the LHC, as well as the advanced status of the pertinent theoretical calculations. 
With the current statistics at the hadron colliders, the workhorse decay channel, at all relevant H 
masses, is H —> WW, followed by W — >> iv, £ — e or \i. Using phase-space singularity techniques, 
we construct and study a plethora of "singularity variables" meant to facilitate the difficult tasks 
of separating signal and backgrounds and of measuring the mass of a putative signal. The simplest 
singularity variables are not invariant under boosts along the pp or pp axes and the simulation of their 
distributions requires a good understanding of parton distribution functions, perhaps not a serious 
shortcoming during the boson hunting season. The derivation of longitudinally boost-invariant 
variables, which are functions of the four charged-lepton observables that share this invariance, is 
quite elaborate. But their use is simple and they are, in a kinematical sense, optimal. 
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It is nice to know that the computer 
understands the problem. But I 
would like to understand it too. 

Eugene Wigner 



I. INTRODUCTION 

Recent data from the LHC [1 on a putative standard 
Higgs boson exclude, at a 95% confidence level (CL), the 
mass domain 127 GeV < M H < 600 GeV (CMS) and, 
with some narrow gaps, 131 GeV <M H < 453 GeV (AT- 
LAS). These results are obtained with full use of the stan- 
dard theory, including radiative corrections which some- 
times constitute the dominant effect. The amplitude for 
Higgs boson production, for example, is largely domi- 
nated by gluon fusion via a t-quark loop and so is the 
amplitude for H — >• 77 decay. 

In the "quantum- level" setting we recounted, it would 
be inconsistent not to analize the LHC data in conjunc- 
tion with the constraints on Mh which follow from the 
profusion of high precision measurements that test the 
standard theory beyond tree level. These constraints 
(and the direct searches [2] at the Tevatron) result in 
M H < 161 (156) GeV at a CL of 95% [3], while the direct 
LEP limit is M H > 115 GeV. 

In mass intervals akin to the one implied by the quoted 
constraints CMS finds a 1.9cr excess of events -that could 
be an indication of a Higgs signal- at Mh = 124 GeV and 
ATLAS a 2.5a one at M H = 126 GeV [TJ. In the current 
broad mass range(s) of the searches, the corresponding 
"local" significances are somewhat larger pQ , but have no 



rigorous statistical interpretation. 

For a standard Higgs boson of mass Mh > 140 GeV 
the branching ratio for the decay H — » WW is the dom- 
inant one. Below this mass and above the LEP limit, 
the winner is H — >• 66, a process beset by terrifying 
backgrounds at a hadron collider. The branching ra- 
tio W — » qq is one order of magnitude larger than the 
one for W —> lis, £ = e, \i. But a light-lepton signal 
is much "cleaner" than that of a quark- generated jet. 
This makes the chain H -> W+W~ , W ± H> t±v the all- 
mass workhorse at a hadron collider. In brief, we refer 
to this process as H — > WW, including the "off-shell" 
M H < 2 M w case, often dubbed H H> WW*. 

The obvious problem with the H — > WW channel is 
that Mh cannot be reconstructed event by event, as a 
lot of information escapes detection with the unobserved 
neutrinos and, at a hadron collider, also with the un- 
observed hadrons that exit "longitudinally" close to the 
beam pipe(s). This makes taming the workhorse almost 
an art, not only a science. The formal and theoretically 
optimal singularity variable procedure to deal with this 
kind of incomplete information is summarized in [4j [5] 
and exploited for the hadron-collider production of a sin- 
gle W in |5]. We shall see that, for the H — >> WW process, 
the situation is much more challenging, mainly because 
two missing neutrinos are many more than one. 

The other crucial obstacles in the process we study are 
the large backgrounds with kinematics akin to those of 
the signal. The main and irreducible one is the direct 
non-resonant production of W pairs by qq annihilation. 
The next most relevant one is ti production, which also 
results in W pairs. For simplicity we shall illustrate our 
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theoretical results only for the chain W —> ev, W —> jiv, 
for which the "Drell-Yan" background is not a problem. 

The analysis tools used to deal with the H — » WW 
channel range from a simple "cut and count" approach 
to "matrix-element" techniques and avant-garde neural 
networks or "boosted decision trees" [1 . There is no 
question that in the long range the methods that input 
and utilize the largest amount of information are likely to 
be the most powerful ones. Whether this is also the case 
at an exploratory "Higgs-hunting" stage is more doubt- 
ful. Here we shall explore a "copy and paste" avenue of 
intermediate sophistication: the derivation of singular- 
ity variables -functions of the observable momenta and 
of Mjj- whose measured histograms are to be compared 
(in one or more dimensions) with pre-prepared templates. 

For the production of a Higgs boson at a hadron col- 
lider, followed by the decays H -» W + W~ , W ± -» t^v, 
we shall limit our discussion to the distributions of vari- 
ous functions of the charged lepton three-momenta k and 
I. The treatment of the 2D transverse momentum of the 
final state hadrons, pr, deserves a separate paragraph, 
the next. The use of other observables, such as the num- 
ber of jets, is beyond our scope. 

Two practical problems are that k and I are measured 
to much higher precision than pr and that the formulae 
for the H — >• WW singularity variables are much more 
complex for pr ^ than for p T = 0. We deal with both 
problems by setting p T = in our theoretical expressions. 
This is less cavalier than it seems. The transverse mo- 
mentum of a Higgs boson in a given event is c± —pr- For 
a given ansatz Mh value, its observed lepton momenta 
can be Lorentz boosted closer to the p T = frame. The 
precise boost would require knowledge of the boson's lon- 
gitudinal momentum, ps. But, typically, p\ <C 2M^, it 
is a fair approximation to neglect ps- More importantly, 
the singularity variables for p T = are very useful even 
in the analysis of events with pr ^ 0, even if one does 
not boost the events back closer to the pr = frame, and 
even if one is also dealing with the quoted backgrounds, 
whose W pairs do not have a fixed invariant mass. 

Let the lepton momenta be k = {kr,ks} and I = 
{lr,h}- Because of the rotational symmetry along the 
beams' axis, the six-dimensional observable space {fc, 1} 
is in practice just five-dimensional. One possible choice of 
variables is the set (k + l) 2 , the invariant mass of the lep- 
ton pair; kr, It-, the moduli of the transverse momenta; 
and kr • It, or the familiar Aip — arccos[/cT • W/ikr It)]- 
All four of these "transverse" observables are invariant 
under longitudinal boosts along the beams' axis. The 
remaining variable, for instance k$ + Z3, is not. 

We shall derive two types of singularity variables: 
those which do -or do not- depend only on transverse 
observables. Transverse variables are preferable, in that 
they are insensitive to the significant uncertainties associ- 
ated with the (longitudinal) parton distribution functions 
(pdfs). In practice the uncertainties are to a modest ex- 
tent reintroduced via the angular coverage limitations of 



an actual experiment, which are not invariant under lon- 
gitudinal boosts. The histograms of singularity variables 
that are not longitudinally invariant do depend on the 
pdfs, but, particularly during a Higgs-hunting or initia- 
tory epoch, this is not a serious limitation. 

In the problem at hand, the quintessential function 
of transverse variables -in the sense of its ability to tell 
apart signal from backgrounds- is Aip, the angle between 
the charged leptons in the transverse plane. In Fig. [I] 
for comparison with our coming results, we recall this 
fact by showing the (arbitrarily normalized) shapes of 
signal and background distributions for two examples, 
with M H = 500 and 120 GeV. As is well known, the V- 
A nature of the weak current and the specific spin zero 
nature of the W pair in signal events, favours collinear vs. 
anticollinear leptons. The effect weakens as Mjj increases 
and the leptons are boosted away from each other. 

The simulations of Fig. [T] as well as all others in this 
note, were made with use of the PYTHIA6 event gen- 
erator [6 . They are for the H — > WW, W — > ev, 
W — >• \iv channel, with leptons of transverse momentum 
greater than 15 GeV, and satisfying the pseudorapidity 
cuts 77(e) < 2.5, r)(ii) < 2.1. 



as 



0.5 1.0 1.5 2.0 2.5 3.0 



0.5 1.0 1.5 2.0 2.5 3.0 



3500 
3000 
2500 
2000 
1500 
1000 
500 



0.5 1.0 1.5 2.0 2.5 3.0 



0.5 1.0 1.5 2.0 2.5 3.0 



Figure 1: Distributions of Acp. Top row: M H = 500 GeV. 
Bottom row: Mh — 120 GeV. Left column: Comparison be- 
tween the shape of the (yellow) signal distribution and that 
of the (blue) WW background. Right column: Comparison 
of signal with the (orange) ti background. 



Our goal is two- fold. First and foremost, to derive the 
complete set of phase-space singularity variables (func- 
tions analogous to Aip) for the process at hand. Second, 
to illustrate with examples their potential phenomenolog- 
ical usage. At least at low Mh values, Acp is more heavily 
dependent on dynamics than on kinematics. The singu- 
larity variables we shall derive are the other way around. 
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Individually, several of them are nearly "as good" as Acp 
in disentangling a signal from the backgrounds. The en- 
semble of their distributions, particularly in conjunction 
with Aip itself, should be a powerful and relatively simple 
tool to search for a Higgs boson, which, in the sense of 
signal kinematics, is guaranteed to be optimal. 

Since we investigate a plethora of singularity variables, 
this paper is long and detailed. The reader mainly inter- 
ested in results may be well advised to start reading it 
from the end: 9XII and 9XII 



II. OUTLINE 

The simple example of single-W production is used in 
§III| to clarify what singularity conditions and singular- 
ity variables are. After posing the formal problem in §IV| 
we proceed in ^V|to solve it in the center of mass (CM) 
reference system of the Higgs boson. There are two rea- 
sons for this. First, it is a necessary intermediate step 
in the theoretical derivation, in §VI[ of the general case 
with a boson which is not at rest. Second, the "approxi- 
mation" of a heavy particle made at rest by gluon-gluon 
or qq fusion in a hadron collider is not so bad, since the 
quark and gluon pdfs are fast falling functions of their 
fractional momenta. Our results will reflect this fact. 

We shall have to deal with the case Mh < 2Mw, in 
which at least one of the Ws is off-shell; the relative 
probability for both of them having an invariant mass 
significantly different from Mw is small. In § VII| we ex- 



plain the simple way in which we treat this case. Further 
details of our data analysis not already discussed in the 
introduction are given in §VIII| 



We analize MC-generated data in [IX in the theoretical 
approximation of a Higgs boson made at rest. To some 
extent, this section is a "warm-up" for the general results 
wherein we lift this approximation, discussed in ^Xj to 
which the reader interested in the most powerful results 
may prefer to jump. 



A summary of results is given in [XI Our data anal- 
ysis is not as thorough as the theoretical one, it is only 
meant to illustrate our points. But it suffices to reach our 
conclusions, which, naturally, are drawn in the last sec- 
tion. A very formal but important step in our theoretical 
analysis is relegated to the Appendix. 



III. SIMPLE SINGULARITY VARIABLES 

Our main result is the theoretical derivation of the 
phase space singularity variables and singularity condi- 
tions for the process H ->> W + W~ , W ± -» i^v. To 
understand these concepts it is easiest to recall a simpler 
problem: the analogous one for single- W production at a 
hadron collider, followed by the same leptonic decay. In 



this case, the singularity condition [5 is Yi T — 0, with: 

z t (mJ t ,Pt) = 

M 4 -4M 2 (f T -p T +^)+4 [(l T -V T ?-l 2 T V 2 T ] (1) 
Of the four M-roots of Yi T = 0, one is not unphysical: 



M T (/ T ,p T ) = + y2 [|Z r |b + Z| T +Z T -(Z T +p T )J, (2) 

which reduces to Mt = 2 \l T \ for p T =0. The function in 
Eq. Q is the habitual M 2 originally derived in [7J |8] . 

The result of Eq. is obtained by projecting the full 
phase space (which includes the neutrino momentum) 
onto the observable phase space. The function S T is 
a singularity variable which -for a general non-singular 
event- is a measure of its distance to the nearest singular- 
ity at the singular St = border of the projected space. 
In Eq. the mass of the W appears in two ways: the 
physical Mh is imprinted in the data and also appears as 
an implicit "trial" mass M — >• M. in the equation. In the 
Mt singularity variable of Eq. ([2| Mh is only reflected 
in the observables. In applying the phase- space singular- 
ity approach to our two-W problem, we shall encounter 
both types of singularity variables. 

In the single-W case one can refine the result in the 
sense of finding the optimal singularity variable, that 
which would result in the most precise measurement of 
Mw [5 . In the two-W case this is not worthwhile, as 
there are decay channels, such as H — >> 77, H — » ZZ; 
Z — >• for which the mass is reconstructible. 



IV. THE FORMAL PROBLEM 

Back to the H -> W+W~ , W ± -> l±v process, let y 
and x, respectively, be the four-momenta of the neutrinos 
accompanying the charged leptons of four-momentum k 
and /. The full information relevant to the reconstruction 
of the boson's mass for a signal event is embedded in the 
kinematical equations: 



Ei 
E 2 
E 3 
Ei 



x A = 

y 2 = 

21 -X = 

2k-y = 



w 



w 



E h ^2{l + x)-{k + y)=M%- 
E e ^ fci + y\ + h + xi + pi = 
E 7 ^ fc 2 + 2/2 + h + x 2 + pi = 



2M^ 



(3) 



where we have made the approximation I 2 = k 2 = for 
the charged leptons and -fleetingly in error for the VKVK* 
case- set the masses of the two Ws to their central val- 
ues. There are 9 unknowns (2 neutrino four-momenta 
and Mh) and only 7 equations. In spite of this, is there 
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a systematic way to extract the kinematically most strin- 
gent information on M#? This is the problem to face. 

Consider the 14D space of the components l,k of 
the three-momenta of the two (approximately massless) 
charged leptons and the four- moment a y of the two 
neutrinos. For a fixed M#, the seven Equations (|3| de- 
fine a 14 — 7 = 7D manifold, the phase space. This surface 
is to be projected onto the 6D hyper-plane of observable 
three-momenta. The points in the full phase space that 
project onto the boundary of the 6D space of observables 
are singular: at such points one or more of the invisible 
directions are contained in the tangent plane to the full 
phase space [H [5] , and a tangent to a surface is singular 
in that it "touches it" at more than one single point. 

The equation, E(Z,fc, M^;M^) = 0, describing the 
boundary of the projected phase space is a singularity 
condition. A general event (i.e. specific values of I and k) 
is non-singular and its corresponding value of E is, once 
again, a measure of distance to the E = singularity. 
The shape of the distribution of the values of the singu- 
larity variable E is sensitive to the unknown mass Mh in 
a manner that allows one to extract its true value, be it 
physical or Monte Carlo (MC) generated. 

The formal modus operandi to obtain singularity vari- 
ables is summarized in [4 and discussed in detail in [5]. 
We recalled that at a singularity one or more of the invis- 
ible directions are contained in the tangent plane to the 
full phase space. The general condition for this to happen 
is that, in the space {z} = {x,y} of invisible directions, 
the row vectors of the Jacobian matrix = dEi/dzj 
(with the row index i running along the number of equa- 
tions and the column index j over the number of invisible 
coordinates) be linearly dependent. In other words, at a 
singularity, the rank of Jij must be smaller than its rank 
at nonsingular points. 

There are 7 equations and 8 invisible directions in 
Eqs. pj. The vanishing of the Jacobian (a 7 x 8 
matrix) entails 8 conditions: the nullification of all 7 x 7 
minors. Two of these minors coincide, up to their sign, 
with two others. Moreover the sums of two pairs of mi- 
nors are of the forms DSq, DSs, with 



D = det(Z, x, /c, y) 

So = k + Zo + + ?/o, 

S3 = ks + h + x 3 + y 3 . 

Given that <So > 0, one condition is: 

det(Z, x, /c, y) = 0, 



(4) 



(5) 



that is, the coplanarity of the four lepton four-momenta, 
equivalent to 



y = a I + j3 x + (7 — 1) k. 



(6) 



Introducing this into the 8 original minors, it is easy to 
see that they all vanish provided that 
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(a - 0) (Z 3 x - lo x 3 ) 



a (k 3 l - k h) + P (A3 ^0 - k #3) 



(7) 



The transposed Jacobian matrix, with use of E6 and 
E7 of Eqs. (|]), is 



J(5 ,5 3 ) 



/ 


xo 





lo 





So 


0\ 




-Xi 





-h 








1 




-x 2 





-h 








1 




-x 3 





-h 





-S3 










yo 





k 


So 










-yi 





-h 





1 







-2/2 





-k 2 





1 


V 





-ys 





-k 3 


-S3 


0/ 



, (8) 



where the functional dependence of J on <So,3 has been 
made explicit for later convenience. 

It turns out to be very useful to study the behaviour of 
the 7x7 minors of J under longitudinal Lorentz trans- 
formations, the boosts along the axis "3" of the proton 
beams. To proceed recall that, for the reasons stated 
in the Introduction, we are setting the hadron momenta 
Pi = P2 = 0. Next, parametrize an event in the usual 
Cabibbo-Maximovich manner [9], illustrated in our no- 
tation in Fig. [2] That is, consider the lepton momenta 
as if both W bosons were at rest, boost them by the 
antiparallel motion of the Ws in the H rest system and 
finally boost the Higgs boson longitudinally left or right 
along the beams' axis: 



I = L{y H , ft p ) L(y, ft) (M w /2){1, ft J, 
x L(y H , n p ) L(y, ft) (M w /2){1, -nj, 
k L(y H , rip) L(y, -ft) (M w /2){1, n k }, 
y L(y H , Hp) L(y, -ft) (M w /2){1, -n fe }, 



(9) 



where n&^n^n and ft p are unit vectors, L(y,n) is a 
Lorentz boost along ft with velocity f3 = tanh(?/), y = 
arccosh(7), 7 = Mh/(2 Mw), and analogously for the 
longitudinal boost along n p , of rapidity yn- 
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Figure 2: Top: The H WW, W ->> iv process in a 
Cabibbo-Maximovich parametrization 9 . The vectors ni and 
fik are the directions of the charged leptons (a /i + and an e~ 
in this illustration) in the respective rest systems of their par- 
ent Ws. The overall WW system, shown here at rest, is to 
be boosted along the direction n p of the gluon or qq pair that 
fuse to produce the W pair, resonantly (for the H signal) or 
not (for the irreducible background). 



Label mj, j = 1 to 8, the 7x7 minors of J in Eq. (|8|, 
lacking the row 9 — j of J. Under a longitudinal boost 
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L(yH:^s): they transform as rrij \-> rrij, with: 

mi = j H [mi + /3 H m s + /3h <$ 3 £>], 

fhi = rrii, i = 2, 3, 6, 7, 

m 4 = Jh [m 4 + p H m b - f3 H S D], (10) 

m 5 7# [m 5 + m 4 + /3# S 3 

m 8 = [m 8 + /3h mi - /3# <S £>], 

where we used the definitions in Eqs. Q. 

The conditions rrij = 0, Vj imply that D = 0, and 
consequently that m^ = 0, Vj. Thus, we reach a cru- 
cial point: the general singular configurations can be ob- 
tained by boosts of the ones in the boson's rest system. 
This is one of the reasons why we pause to study this 
latter simpler case. 



V. LESSONS FROM A GLUON COLLIDER 

A standard Higgs can be made in various ways, with 
top-mediated gluon fusion being the dominant mecha- 
nism up to very high M# values. The gluonic pdfs, as 
well as those of the other partons, are fast-falling func- 
tions of their fractional momentum. This implies that 
Higgs bosons are made with a narrow distribution of ra- 
pidities, centered at yn = 0. The same is true for the 
backgrounds to the H — » WW channel, e.g. non-resonant 
pairs of relatively heavy objects, such as VK-bosons, are 
also made with a moderate collective motion. Thus, the 
approximation of a monochromatic "gluon collider" (or 
qq collider) is a good starting point for our analysis. 



A. Derivation of the singularity conditions 

In the WW center of mass (CM) system an extra work- 
ing condition is to be added to Eqs. (J3|: 



Ccm ^ S 3 = h + x 3 + k 3 + y 3 = 0. 



(11) 



and the Jacobian is now J(<So,0), with J as in Eq. 
Since So > and the fifth column of J(So,0) is propor- 
tional to So, it suffices to consider the vanishing of the 
eight 7x7 minors of J(l, 0), of which only four, e.g. m^, 
i = 1, ...4, are independent modulo D. 

In the CM system, let E = M H /2 denote the W's 
energy and P the corresponding momentum modulus. 
The four-momenta of the individual Ws in the notation 
of Eq. ([9} are: 



p wl = {E,+Pn} 
P wk = {E,-Pn} 



(12) 



and it is convenient to put the neutrino's momenta in the 
form x = p wl -l,y = p wk - k. 



The conditions x 2 = y 2 = now read 



2El ~- 
2Eko 



Mi 



w 



Mi 



w 



2Pl-n 
-2Pk-n. 



(13) 



Stepping back to Eq. ([8| and introducing the explicit 
lepton four-momenta in the minors of J(1,0) and in 
det(Z, x, fc, 2/), the vanishing of the results requires, in par- 
ticular, that det(Z,fc, n) — 0, that is, the 3D coplanarity 
of Z, k and n and, consequently, of all four lepton three- 
momenta. We may write 



n = (a/lo)l + (b/ko)k. 



(14) 



Gathering results and imposing n • ft = 1, one may 
express a, 6, x and y as functions of Z and k. Two families 
of CM critical configurations are obtained. They differ 
by the sign of 5 in 

S 2 = (4k 2 - 2M H k + M^)(4Z§ - 2M^Z + M 2 ,) (15) 



and satisfy: 

k-l = 

Zo[ 



-4,E 2 k lo + 2E(k 



Z )M^ - 1V 



Mi 



a 



4 p2 

2Ek (k ■ 1 + k l ) + {kl + k- l)M$y] 



2[{k-lf-klll]P 
k [2El (k ■ 1 + k l ) ~{ll + k- f)M&] 



2[{k-lf-klll]P 



(16) 



Substituting these expressions into mj, i = 1, ...4 one 
finds that vanishes automatically. The others inde- 
pendent minors acquire the form 



mi = (k 2 h-k 1 l 2 )Af/V, 
m 2 = (ksh-k^Af/V, 
m 3 = (k 3 l 2 -k 2 h)M/V, 



(17) 



where M and V are lengthy functions of k and Z. 

There are two alternative ways to satisfy rrii = Vz. 



One of them is to let all three parenthesis in Eqs. (17) 



vanish simultaneously, tantamount to imposing k oc Z, a 
specific case of the condition to be obtained anon from 
the second alternative: J\f = 0. Eliminating the sign 
ambiguity of S yields a first requirement for an event to 
be singular, C = 0, with 



k -2M^°M^ 
C\ = k() — Zo, 

C 2 = 2k l M H - M^(ko + Z ), 
C 3 = 4M^(fc Z 3 + hlo) [2M H (kok 3 + Z Z 3 ) + 
2(k - ZoXMs - k 3 l ) - M 2 H (k 3 + Z 3 )] + 
M^(k 3 + l 3 ) [M H (k 3 + Z 3 ) - 2(k - lo)(h - Z 3 )] 
4M H (2k - M H )(2l - M H )(k l 3 + k 3 lo) 2 



lH ,, w (M^-4M^) 3 . 



(18) 



The non-trivial vanishing of C 3 implicitly presupposes 
det(Z,/c,n) = 0. Up to non-vanishing overall factors, a 
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second requirement for an event to be singular is this 
coplanarity condition, squared such as to eliminate the 
sign of 5: Co = 0, with 



Co = 2M^(hk - 2k l ) [21 -k- M H (k + Z )] (19) 
- M H (hk - 2k lo) 2 + [21- k - (k + lo) 2 } , 
where l-k has its customary Minkowskian meaning. 

B. Questions of nomenclature 



For a singular event the values of C in Eq. (18) and 
Co in Eq. (19) must both vanish. Given the form of 



C, there are three nontrivial ways for this to happen: 
Ci = Co = 0, i = 1 to 3, which we call complete singularity 
conditions. Of these, only Cs = Co = guarantees that all 
minors of the Jacobian vanish. The other two conditions, 
Ci = Cq = 0, i = 1 to 2, are mock singularity conditions, in 
a sense occasionally used in mathematics, that is, they 
do not satisfy all wanted conditions, but are useful for 
one's purposes. As it turns out, even the four partial 
singularity conditions Cj =0, j = to 3, are of interest. 

We choose Co as the example to make our next linguis- 
tic points. Consider a real or MC-generated event due 
to the production and decay of a Higgs boson. Its corre- 
sponding value of the Co function in Eq. ( p~9| ) -a (partial) 
measure of distance to the Co = singularity- depends 
on the Higgs boson mass in two distinct senses. The 
first is that ko, lo and /care contingent on this "input" 
mass. The second is the explicit Mh in the expression 
of Co, which is a variable that one may -naturally- vary 
at will. To emphasize this point, we label this analyst's 
mass calligraphically: Mh —> M. 

It is convenient to rescale and rewrite Co as: 

Co = -M|S (A^)/4 
Z (M) = (M - M+)(M - M_) 
2M W 



M± 



Ml 



M w (ko + lo)±M M \ M } 



r 2 



Mh E = 2(k l Tkl ), 



(20) 



where M± are the non-zero roots of Co — and M^ E are 
the Minkowski and "Euclidean" masses of the (approx- 
imately massless) charged lepton pair. Notice that So 
depends on the implicit variable M, while its roots, M± 
do not. That is why we refer to So and M± with different 
symbols, even though their distributions are in all cases 
diagnostics of the value of the real or simulated Higgs 
boson mass (we reserve the nomenclature "M" for all 
singularity variables of the later kind). Implicit masses 
become theoretically inevitable in cases for which, unlike 
for So, the M. roots cannot be made explicit. 

Functions of an implicit mass M., such as Eq, are also 
singularity variables. They vanish at singular points of 
phase space, iff the correct choice M. = Mh has been 
made, with Mh the physical or Monte Carlo "truth" . 



Partial and complete singularity conditions and 
variables 



The singularity condition C = of Eq. (|18| can be 
satisfied in various ways. Two of them (Mh = and 
Mh = 2 Myy) are of little practical relevance. Two others 
correspond to the naive-looking observables 



Mi 
M 2 



\ko ~ k\ 



2 ^0 + ^0 
w 



2ko lo 



(21) 
(22) 



The remaining possibility is C3 = in Eq. ( |18| ), a 
cubic polynomial in the Higgs boson mass. In analogy 
with Eq. (20) for Co, we introduce its roots: 



C 3 = F S 3 (7W) 
T, 3 (M) = (M- M X ){M 
F = 4(k 3 l + k 1 3 ) 2 



M 2 )(M-M 3 ) 



(23) 



where the explicit forms of Mi are lengthy. 

It is not useless to rewrite Eqs.( 2Q|22|23 ) in the form: 



S ee Sign(S ) |£ | 1/2 

El EE Ml 

E 2 = M-M 2 

E 3 ee Sign(E 3 ) |S 3 | 1/3 



(24) 



This is because to construct true singularity variables 
that reflect a complete set of singularity conditions we 
must introduce a measure of the distance between a data 
point (given values of k and I ) and one of the three 
center-of-mass singular manifolds: the points {0,0} of 
the planes {C j = 0} n {C = 0}, i = 1 to 3. With the 
help of Eqs. ( |20|24 ) we define the following quantities 
with unit mass dimension: 



D 2 = (e^ + 5 



D* 



(Eo + E|) 



1/2 
1/2 
1/2 



(25) 



The functions Di(A4) are the full set of complete center- 
of-mass singularity variables for the case at hand. 



D. From Algebra to Geometry 

An advantage of the approximation in which Higgs 
bosons would be produced at rest is that the locus of 
the singular points in the observable {/c, /} phase space 
can be visualized. Let c^ = cosA(/? = k T • It/(^tIt)- 
The singular phase space is shown in Fig. [3] in the vari- 
ables {&o,/o? c ^}? m an example wherein we have chosen 
Mh = 2.5 in Mw = 1 units. The closed surface in the 
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three subfigures is the coplanarity condition Co = 0, see 
Eq. (19). The thick lines in the top and middle figure 
correspond to the singularity conditions Co — C\ = 
and Co = C2 = 0, see Eqs (18). The last figure partly 
describes the singular phase space Co = C3 = for the 
choice /C3 + /3 = 0; the complete space would be the direct 
product of this latter line in {^3,^3} space with the thick 
line in the figure. 

The Co = Cs = singularity condition, as one varies 
&3 + ^3, covers all of the Co = coloured surface of Fig. [3] 
This reflects the fact that the other two conditions are 
mock, and of zero measure relative to Co = Cs — 0. 



VI. BACK TO A HADRON COLLIDER 

The derivation of a singularity variable for the more 
realistic case of an H boson of rapidity yu 7^ is akin to 
that of the yn = case, requiring only one extra step. 
Naturally, this is to start by applying the Lorentz boost 
L(yH,n p ) to the W momenta of Eq. (12), to obtain: 



Pwi = i cE + s P n 3, +Pni, J rPn 2 , cPn 3 + sE} 
Pwk={ cE ~ ^Pn 3 , -Pn u -Pn 2 , cPn 3 - sE} (26) 

where c = cosh(^) and s = sinh (?/#). Following pre- 
cisely the same steps as in j VA[ one concludes that the 
partial singularity conditions are C[ = 0, i = l to 4, with 



c" 2 



C'o 



2£ Ci(l' , l 3 , k' , k 3 ) 

k + k 3 - l - l 3 + (k - k 3 -l + h) £ 2 

2eC 2 (l' ,l' 3 ,k' ,k' 3 ) 

(k + k 3 ) (l + k) M H - ... 

+ (k -k 3 ) (l -h)M H ^ 

Aec 3 (i' ,i' 3 ,k' ,k 3 ), 

4(ko + k 3 ) 3 (l + h) 3 M H - ... 
+4(k -k 3 ) 3 (l -l 3 ) 3 M H e 2 
^ 4 C (l' ,l' 3 ,k' ,k' 3 ) 
-(k + k 3 ) 2 (l + h) 2 M 2 H + ... 
-(k -k 3 ) 2 (l -l 3 ) 2 M 2 H e 



where 



(27) 



(28) 



f = c + s = e VH 
l'o = c ^o - 5/3, I3 = els ~ slo, 
k'o = cko — 5/C3, k' ?j — cks — sko 



In C2307 whose expressions in terms of unprimed mo- 
menta are easily obtained and lengthy, we have only given 
the first and last term in their expansion in £, which are 
sufficient to specify their mass dimension and their grade 
as polynomials in £, two numbers that we shall need. 

To obtain longitudinally boost-invariant results anal- 
ogous to the ones in Eqs. (25) one must eliminate the 




(19) 



Figure 3: Visualizing the singular CM phase space, 
closed surface in all figures is Co = 0, see Eq 
thick lines in the top (middle) figure are the singu 
ditions Co = Ci = (Co = C 2 = 0), see Eqs <\18) 



The 
The 
arity con- 
The 



unknown boost parameter £ between the pairs of poly- 
nomials {Cj(£), C*q(£)}, j = l to 3. The first and simplest 



thickest black line in the bottom figure is the singular phase 
space Co = C3 = for the particular choice + £3 = 0. The 
entire illustration is for Mw — 1,Mh = 2.5. The horizontal 
axes are kr and It, the vertical one is = kr • It /(kr It)- 
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of these results, for j = 1, is the singularity condition 
Ai = 0, with 



Ai oc M 2 E, 

M = A-Bl- k, 

E = [4A 2 E 4 + B 2 (M^ - 2P 2 Z • fc) 2 
-4A££ 2 (M^ + 2P 2 Z-/c)], 



(29) 



(l T -k T + l- k) 2 



k^ Irp 



B = 2 (Z T • fc T + Z • fc) - k 2 T - 1%, 

where E = M./2 is the energy of a W in the rest system 
of a Higgs boson of trial mass M. and P is the corre- 
sponding momentum. We have followed our convention 
to label M the singularity variables that do not depend 
on M., and E (and now A) those which do. Notice that 



Ai, by construction and demonstration, is a function of 
longitudinally boost-invariant observables. 

The derivation of analytical results for the remaining 
polynomial pairs is not as simple as it was for {Cq,C(}, 
C[ being merely quadratic in £. The expressions for 
C' 2 3 o are polynomials in £ of degrees 4, 12 and 8, re- 
spectively. The condition for two polynomials 5^=0 
and Y^j=o & to vanish simultaneously (to have com- 
mon roots) is called their resultant, and is a sum of prod- 
ucts of powers of ai and bj. The resultant of C[ = 
and Ch = is the condition Res{C{,Co} = Ai = 0, see 
Eq. ( |29| ). The number of terms of a resultant grows very 
rapidly with ra x n, it is 95 for (ra, n) = (2,8), 4970 for 
Res{C2,Co}, for which (ra, n) — (4,8). For this case, af- 
ter considerable simplifications, the singularity condition 
is A2 = 0, with 



Res{C 2 , C' } = A 2 oc £ 4 k% Z§, (l-k) A (M$(l§E A k 2 T l 2 T + (-4(t T -k T ) 2 + 2t T -k T M 2 v + Z 2 M 2 , + k\ (4Z 2 + M 2 ,)) 2 

- SE 2 (2k A l 2 T + l T -k T l 2 T (-2l T -k T + M 2 ,) + /c 2 (-2(Zt-^ t ) 2 + T T -k T M^ + 2Z 2 (Z 2 . + M 2 ,)))) 
+ 41 -kM$ (-32 E 6 k^ l 2 T + M^(8(Zt-^ t ) 3 + 2l T -k T M^ + Z 2 M 2 ^ 2 + M 2 ,) 

+ 4(4Z 2 + M 2 ,) - 4(Zt-^ t ) 2 (Z 2 . + 2M 2 ,) + 4(-4(Zt^ t ) 2 - 8Zr-Mr + 4/ t + ^M 2 , + M 4 ,)) 

- ^ 2 M 2 7 (-8(Zt-^t) 3 + 8/t-MtM 2 , + 4(Zt-^t) 2 (-5Z 2 + M 2 ,) + /c 4 (20Z 2 + M 2 ,) + Z 2 M 2 ^ + 2M 2 ,) 
+ 2/c 2 (-IO(Zt-^t) 2 + 10Z 4 + 13Z2.M 2 , + M 4 , + U T -k T (l 2 T + M 2 ,))) 

+ 4£ 4 (4A;|4 + l T -k T l 2 T (-Al T -k T + 3M^) + 4(-4(f T -£ T ) 2 + 3l T -k T M^ + 2Z 2 ^(2Z 2 . + 7M^)))) 
+ 16^ 4 (Z-^) 4 (M^ + 4^ 2 P 4 )(M^ + 4Z 2 P 4 ) + 4(Z^) 2 M^(-4Zt^ t M 2 ,(2^ 2 + M 2 ,) 
x (3£ 4 (fe 2 + Z 2 ) - 4£ 2 (& 2 + Z^M 2 , + M^(4 + Z 2 + M 2 ,)) 

+ 4(Zt-£ t ) 2 (4£ 6 (/c 2 + Z£) + + E A M^(-^k 2 T - 8Z 2 + M 2 ,) + 4£ 2 M 4 F (& 2 + Z 2 + M 2 ,)) 

+ M^(M^(Z£ + M 2 ,) 2 - 2E 2 l 2 T M 2 v (l 2 T + 6M 2 ,) + £ 4 Z 2 (Z 2 . + 12M 2 ,)) + 2*£(48£ 8 Z£ + M^(Z^ + M 2 ,) 

- 8£ 6 Z 2 (Z 2 + 15M 2 ,) + £ 4 (16Z 4 M 2 , + 97Z 2 .M 4 F + 6M^) - 2£ 2 (4Z 4 M 4 , + lSZ^M^ + 3M^)) 
+ 4(-16£ 2 Z 2 + M^)P 4 ) - 16£ 2 (Z-*0 3 M^(M# + 3(4 + l 2 T )M^P 2 

+ 6(/c 2 + l 2 T )M^P A + 244Z 2 .M 2 7 P 6 + 32k 2 T l 2 T P* - 2l T -k T M 2 v (E 2 + M^)(M^ + 2(fc 2 + Z 2 )P 4 ))) 



Notice that A2, as was the case for Ai, is a function of 
only longitudinally boost-invariant observables. 

What is the number of terms in ResjC^, Cg}, for which 
the degrees of the polynomials in £ are (ra, n) = (12,8)? 
For ra, n larger than a small integer the resultant soon 
becomes obdurately complex. Not even the number of 
addends in the (monomial) products of its formal coeffi- 
cients is known. Only upper bounds to that number are, 
to our knowledge, published. The tightest one is [TP] : 



S(m,n) = P(ra, n, [mn/2\) 



where, for integer a, 6, c, F satisfies the recurrence 

b 

F(a,b,c) = ^P(a- l,j,c- j), 

with 



F(l,6,c) 



1 if < c < b 
otherwise 



For the polynomial pair {C' 2: Cq}, (ra, n) — (4,8) and 
5(4,8) = 16335 (an overestimate by a factor ~3), while 
for {C^,C£}, (ra,n) = (12,8) and 5(12, 8) = 477174 360. 
This last upper limit is the best current estimate (by 
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expert mathematicians) of the number of terms in the 
expression for the remaining singularity variable we are 
after, in terms of products of powers of the m + n + 2 = 22 
coefficients of £ in the polynomial pairs, each of which is 
a complicated function of Mw, A4, k and /. 

We shall not be discouraged by the mathematical hard- 
ship of constructing explicit algebraic resultants. In anal- 
ogy with Ai in Eq. (29) and given the complexity of A2 



in Eq. (30), we shall simply define: 



A 2 ee Res{C' 2 (0,C' (0} 
A 3 == Res{C 3 (0,^(0} 



(31) 



and find, event by event, the resultant numerically. The 
coefficients of the powers of £ in C' 2 3 being -for a given 
event- numbers as opposed to symbols, this is doable and 
-for the computer- trivial. 



The formal proof that the resultants in Eqs. (31 ) ought 



to be boost-invariant is given in the Appendix. 



VII. DEALING WITH THE M H < 2M W CASE 



(M 2 ) 



2\ 0.75 



Mw 0.7 



0.65 
0.6 
0.55 



100 



M H (GeV) 



150 



Figure 4: The average squared mass of Eq. ( 32 ) , as a function 
of M H , in the process H -> WW*, for M H < 1.93 M W - The 
continuous line is the leading-order calculation in the standard 
model. The dashed line is the approximation of Eq. (34). 



In an H — >> WW* process followed by leptonic decays 
of both Ws, there is no way to assign a mass, M*, to 
the W which is putatively off-shell, even for a fixed M#. 
Moreover, there is no deterministic way to decide which 
W was approximately on-shell. Finally, except close to 
the Mh = 2Mw threshold, the theoretical distribution 
of off-shell masses, dT/dM*, is very wide. To confront 
this situation we choose to analize this case by assigning 
to both Ws an adequately averaged squared mass: 



(M 2 ) = 



(32) 



Since the non-observation of two neutrinos results in wide 
distributions for all observables, there is very little differ- 
ence between using this prescription and other sensible 
ones, such as substituting the average in Eq. (32) by 
its most probable value. 

For Mh > M\y, up to a few W widths, below the 
two-W threshold, and for the leading order standard- 
model matrix element for the H — >> WW* decay, the dis- 
tribution of W* masses, in the excellent approximation 
of neglecting T^/M^, is: 



dv 



2M^{M H + Ml) 



(M 2 



H 



Mir 



dM* 



M w 



(MS- 

2M^ (M 2 H - hMl 



M 2 W Y 
+ (M 2 H 



■Mi 



(33) 



The corresponding (M 2 ) distribution is shown in Fig. Q 
for the relevant range of Mh values. In this range and 
to a good approximation 



{M\M H )) = ^ 



M; 



H 



1.02456 Mi 



w 



(34) 



also shown in the figure as the dashed line. 



VIII. DETAILS OF OUR DATA ANALYSIS 

We have derived singularity variables only for the sig- 
nal process, not for its backgrounds, and we use the sig- 
nal singularity variables to compare the distributions of 
MC-generated signals and backgrounds. 

We present results only for the H — >> WW, W — » eis, 
W \iv channel and its non-resonant WW and tt back- 
grounds, with leptons of transverse momentum greater 
than 15 GeV, and satisfying the pseudorapidity cuts 
77(e) < 2.5, r/(/i) < 2.1 0. 

Given the delicacies of measuring or simulating (at 
a "reconstruction level") the transverse momentum of 
hadrons, p?, we have not boosted each event to the 
approximate frame wherein the putative Higgs boson is 
transversally at rest. Our MC-simulations are for "gener- 
ator level" events and do not have a pr = requirement. 
No doubt this makes our results look somewhat weaker 
than they might otherwise be. 

The ratios of signal to background yields are fast- 
varying functions of Mh- The selections made by ex- 
perimentalists on the way to focus on signal events are 
many and are also mass-dependent. These are reasons 
why we shall limit ourselves to illustrating only the dif- 
ferent shapes (and not the absolute scales) of the signal 
and background histograms of various singular variables. 

In discussing singularity variables such as T>o(M) of 
Eq. (20) or Ai(A^) of Eq. (29), it is informative to do 
it not only for the correct "guess" M. = Mh-, but also 
for incorrect ones. Naturally, the histograms for a fixed 
Mh and various M contain precisely the same statisti- 
cal information. An experimentalist using an observable 
such as or Ai would deal with data (with Mh not 
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known a priori) armed with a plethora of "diagonal" MC 
templates with M = Mh, with which to compare the 
observed distributions. 



IX. DATA ANALYSIS IN THE CM 
APPROXIMATION 

In this section we sketch a numerical analysis of the 
partial and complete "Higgs at rest" singularity variables 
derived in ^V) Recall that these theoretically-obtained 
expressions ignore both the longitudinal and transverse 
momentum of the Higgs-boson signal to be analized. 



A. Partial singularity conditions 

We start this part of the discussion with the singularity 
variable So, a measure of distance of an event to the 
partial singularity condition of coplanarity: = 0. We 
chose Mh — 500 and 120 GeV as examples of the "true" 
mass of the events in this first illustration. 

The distribution of values of S for 20000 signal events 
generated with Mh = 500 GeV is shown in the left panels 
of Fig. [5] The top left panel is for the correct assumption 
M. — Mh, the two other left panels show comparisons 
with the incorrect assumptions M = 4Mh/5 (middle) 
and M = bMn/^ (bottom). The right panels of Fig. [5] 
show results for Mh = 140 GeV. The top panel is for 
the correct assumption M = Mh- The middle panel is 
for M = 120 GeV and the lower one for M = 160 GeV. 
At Mh = 500 GeV the distribution of is very sensi- 
tive to the boson's mass, as exemplified in Fig. [5] by the 
sensitivity to M. At M# = 140 GeV this is less so. 

The ability of the E distribution to sieve apart signal 
and background shapes is illustrated in Fig. [6j Its left 
(right) columns are for M H = 500 (120) GeV, both with 
A4 set to its corresponding correct value. The top (bot- 
tom) lines refer to the WW and ti backgrounds. In all 
cases we have simulated equally many signal and back- 
ground events, so that the figure reflects the shape of the 
distributions, not their relative weights. At Mh = 120 
GeV the shape of signal and backgrounds are very differ- 
ent, while at Mh = 500 GeV this is less so. 

The conclusions on the ability to distinguish signal and 
backgrounds or different Higgs masses are, as we saw, 
very mass dependent. The rest of the questions to be 
discussed in this chapter are quite insensitive to Mh- We 
shall study them only for the Mh = 120 GeV example. 



The quantity So of Eq. (20) is real, but its roots, 
M±, need not be. For input MC data corresponding to 
Mh = 120 GeV, about 14% of the roots are a real pair, 
the rest being two conjugate complex numbers. The con- 
clusion that the complex roots are useless would be most 
premature. A feature of these roots to be studied ab 
initio is the correlation between their absolute value and 
phase. This is done for the Mh = 120 GeV signal and the 
WW and ti backgrounds in Fig. [7[ where the mass axis 





-600000 -400000 -200000 200000 
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Figure 5: The distribution of H values. The horizontal scales 
are in units of GeV 2 . The left column in for M H = 500 GeV. 
Its top figure is for M — Mh- In its middle (lower) ones the 
result is compared with that of the incorrect M = (4/5) Mh 
(M = (5/4) M H )- The right panels are for M H = 140 GeV. 
The top one is for the correct M = Mh- The middle panel 
is for M = 120 GeV, the lower one for M = 160 GeV. In the 
figures the correct-guess histogram is yellow. 



is the absolute value of the roots (shown once for each 
complex root and for its two values for each real pair). 
The cp axis is the phase of the roots having a non-negative 
imaginary part. We see that the {cp : \M\} correlation is 
weak and the distributions are significantly different for 
signal and background. 

The distribution of absolute values and phases of the 
roots of So, that is the projections of the results of Fig. [3] 
onto the \M\ and (p axis, are shown in Fig. [8] The results 
for signal and WW and ti backgrounds are significantly 
different. Notice in particular how the signal has a much 
higher fraction than the background of events with \M\ 
and (p close to zero. 

The variables M lj2 of Eqs. ( |21|22| ) are akin to M± in 
that they do not refer to an ansatz mass A4. In spite 
of their naivete, these observables, particularly M2, are 
quite good at telling signal from backgrounds. Their 
shapes for an Mh = 120 GeV signal and the WW and 
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Figure 6: Shapes of the So distributions for the signal -the 
yellow tallest histogram- and the WW and tt backgrounds 
(upper and lower rows). The horizontal scales are in GeV 2 
units. The left (right) column is for M H = 500 (120) GeV. 



tt backgrounds are shown in Fig. [9] 

Because the variable £3 of Eq. ( [23] ) has mass dimension 
3, it is convenient to plot its sign-recalling cubic root of 
Eq. (24). This we do in the left column of Fig. 10 for an 
Mh = 120 GeV signal and the WW and tt backgrounds. 
The signal and the illustrated backgrounds are seen to 
result in distributions with similar looks but significantly 
different details. In the right column of Fig. [l0]we show 
the three roots of the cubic equation Yi^M.) = 0, see 
Eq. (23). The taller (yellow) histograms are the Mh = 
120 GeV signal, they are compared with those of the WW 
background (the tt distributions, not shown, differ a bit 
more than the WW ones from the signal distributions). 
The three roots Mi of S3, unlike the roots M± of So, are 
not so useful in telling signal from backgrounds. 



B. Correlations between partial singularity 
variables 



are cor- 
tgrounds. 



A question of practical interest is the extent to which 
the Co and Ci^,3 distributions of Eqs. ( 20|18 
related, for a putative signal, and for the bac 
It can be answered, pictorially, by contemplating 2D his- 
tograms in the three {Co,Ci} planes. In the {Co^Cs} 
case, for which the functions depend on A4 : we c hoose t o 
plot the results in the {So, ^3} plane, see Eqs. ( 20|23 ). 
The singularity, for the correct assignment M. = Mh, 
is at the origin of the plane. For this to be the case 




a *jl 




M(GeV) 



Figure 7: Correlation between absolute value M and the 
phase <p of the roots of So- Top: Mh = 120 GeV signal. 
Middle: WW background. Bottom: tt background. The real 
roots gather along the ip — axis. 
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Figure 8: Shapes of the distributions of \M\ (left column) 
and (f (right column) of the roots of Eo- The top row is the 
M H = 120 GeV signal. The WW and ti backgrounds are 
shown in the middle and lower rows. \M\ is in GeV units. 




Figure 9: Comparison of the shapes of the distributions of M2 
(left column) and Mi (right column) for an Mh = 120 GeV 
signal (yellow) and the WW (top row) and ti (lower row) 
backgrounds. See Eqs. ( 22|21 ). Abscissae in GeV units. 



are shown in Fig. [l2j The three choices appear to be 
comparably efficient at telling signal from backgrounds. 

The {Di,Dj} correlations are illustrated in Fig. 13 for 
a signal with M H = 120 GeV and for the WW back- 
ground, all with M. — Mh. Only the {D2 ,^3} corre- 
lation is strong. In all cases the signal and background 
results are fairly distinct. This is even more so for the ti 
background, results for which we do not show. 



in the two other pairs, we plot results for {Eq,Ei} an d 

{E , £ 2 }, see E <1- (If])- A11 this we do in Fi s- @- 

Two conclusions are to be extracted from the quoted 
figure, after noticing that its horizontal scales for signal 
and backgrounds are not always the same. The signal 
variable pairs are quite correlated, with the exception of 
{So, Si}. The WW background is less correlated and 
its distribution is significantly different from that of the 
signal, for all variable pairs. These statements are more 
so for the ti background, which we have not shown. 



C. Complete CM singularity conditions 

To construct true singularity variables that reflect a 
complete set of singularity conditions we must exploit a 
measure of the distance between a data point (its values 
of k and / ) and one of the three center-of-mass singular- 
ities which, for M = Mh, are the points {0,0} of the 
planes {Eq,E^}, i = 1 to 3. These are the quantities 



X. DATA ANALYSIS BEYOND THE CM 
APPROXIMATION 

We have derived three longitudinally boost-invariant 
singularity variables. The first of them, Ai is alge- 
braically simple and factorizable as A 1 ex M 2 E, see 
Eq. (29). The mass dimension of C[ is 1 and its degree 
in £ is 2. The corresponding numbers for Cq are 6 and 
8, see Eqs. (27). The mass dimension of their resultant 
is 1x8 + 6x2 = 20. The mass dimensions of M and E 
are 4 and 12, respectively. It is therefore convenient to 
discuss the results in terms of M 1 / 4 E 1 / 12 and 



Ai 



- A 1/20 



(35) 



where the root is always real, since Ai is always positive. 



Di defined in Eqs. (25). Distributions of these variables 



For the quoted variables we show in Fig. 14 histograms 
comparing the quite distinct shapes of the distributions 
for a signal of a Higgs boson of mass Mh = 120 GeV and 
the WW and ti backgrounds, for M. = Mh in the case 
of E and A 1 . The distributions of M 1 / 4 and E 1 / 12 -and 
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Figure 10: Left: distributions of the variable E3 of Eq. (24) 



top is the M H = 120 GeV signal, middle is the WW back- 
ground, bottom is the it one. Right: distributions of the three 
roots of E 3 (.M) = 0, for the quoted signal (tall and yellow) 
compared with the WW background. Abscissae in GeV units. 




Figure 11: Top left: Histogram of the values of Eo and Ei, for 
a signal with M H = 120 GeV. Top right: for the WW back- 
ground. Middle row: the same as the first row, for {Eo,^}. 
Third row: {E ,E 3 }. All plots are made for M = 120 GeV 
and all horizontal scales (which do not have the same extent 
in all fi gure s) a re i n GeV units. The relevant definitions are 
in Eqs. pol [23l [24). 



consequently Ai- are similar, since the first two variables 



are correlated. The correlations, shown in Fig. (15), are 



not as strong as one might have suspected on the basis of 
Fig. [14] . The correlations are weaker for the signal than 
they are for the background, except in the high-mass tails 
of the background distributions. 

The ability of E and its factors to tell apart diverse 
masses (120 and 140 GeV in the coming instance) is stud- 



ied in Fig. [16] The top left figure, histogramming E 1 / 12 , 
is for a fixed M H = 120 GeV, with M = M H or 140 
GeV. The other tree figures, for M 1 / 4 , E 1 / 12 and A x are 
for M H = 120, 140 GeV with, in all cases, M = M H . 
The function E(A^) is quadratic in M. so that its roots, 
in analogy with M± in Eq. (20) can be made explicit. 



But they are not very efficient at telling apart signal from 
backgrounds, nor at zooming into a value of M#. Thus, 
we do not show results for them. 



We conclude that the singularity variable Ai and its 
factors are strong boost-invariant tools to tell signal from 
backgrounds, and is not very stringent in constraining the 
value of Mh- 

The mass dimension of C 2 is 3 and its degree in £ is 
4. The corresponding numbers for C' are 6 and 8, see 
Eq. (27). Thus, the mass dimension of A2 = ResjC^, Cq} 
48. In analogy with Eqs. (24), it is 



is 3x8 + 4x6 
convenient to define 



A 2 = Sign(A 2 ) |A 2 



1/48 



(36) 



Results for the distributions of this singularity variable 
are presented in Figs. T7J 18 and commented later. 

The mass dimension of C' s is 7 and its degree in £ is 
12. Recall that thee corresponding numbers for C f are 
6 and 8. The mass dimension of A3 = ResjC^Cg} is 
7x8 + 12x6 = 128. In analogy with Eq. pol), it is thus 
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Figure 12: Histograms of the singularity variables Di of 
Eq. ([25J. Top: D x . Middle: D 2 . Bottom: D 3 . Left: compar- 
ison of the signal for M — Mh = 120 GeV (the distribution 
peaking closer to zero) with the WW background. Right: 
comparison of the same signal with the ti background. Ab- 
scissae in GeV units. 



convenient to define 

A 3 = Sign(A 3 ) IA3I 1 / 128 



(37) 



Results for the distributions of this singularity variable 
are presented in Figs.[l7|[l8j The message of these figures 
is that the variables Ai ? 2,3 are very good both at distin- 
guishing signal and background events and A2,3 are very 
good at pinpointing the mass of a putative signal. 

An interesting feature emerges when some of the his- 



tograms in Fig. 18 are remade with higher statistics and 
resolution, concerning the singularity functions A 2 ,3, but 
not Ai. This is shown in Fig. 
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A very clear narrow 
double peak shape appears for A 3 (lower figure), and a 
hint of a similar structure for A 2 (upper figure). 

The peaks in Fig. [19] reflect individual roots in M of 
A2,3(.M). The function A 2 in Eq. (30) is, after elimina- 
tion of the overall E A = (A^/2) 4 factor, still a polynomial 
of sixth degree in M 2 . The A 3 (M) resultant in Eq. (folk 




Figure 13: Histograms of the correlations between the singu- 
larity variables Di of Eq. ( [25] ) . Left: the signal for Mh = 120 
GeV. Right: WW background. Top: {D U D 2 }. Middle: 
{Di,D3}. Bottom: {^2,^3}. Horizontal axes is GeV units. 



has an even more intractable degree in A4 2 : twenty-two. 
Thus, their roots can only be extracted numerically event 
by event, a rather laborious task, which we postpone. 



A. Correlations 

The longitudinally boost invariant singularity variables 
Ai,2,3 have correlations similar to the ones between ^1,2,3 
that we showed in Fig. [13] They are shown in Fig. [20j 
for an Mh = 120 GeV signal. Once again, there are 
significant but not extreme correlations. Moreover the 
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Figure 14: Hist ogra ms of the singularity variables M 1//4 , E 1 / 12 
and Ai of Eq. ([29}. Top: M. Middle: E. Bottom: Ai. Left: 
comparison of the shape of the signal for A4 — Mh = 120 
GeV (the distribution peaking closer to zero) with the shape 
of the WW background. Right: comparison of the same signal 
shape with that of the ti background. Abscissae in GeV units. 



correlated histograms are quite different for the signal 
and WW background. This is even more so for the ti 
background, which we do not show. 



In Fig. 
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we illustrate the correlations between Ai ? 2,3 
and Aip for the signal and the WW background, to which 
the ti background again is in this sense similar. For the 
relatively light Mh = 120 GeV signal shown in the fig- 
ure, as expected, signal and background densely populate 
very different regions of phase space. 



XI. SUMMARY OF RESULTS 

With an eye on potential practical usefulness, let us 
call "good" the singularity variables that do an efficient 
job at focusing on the correct value of Mh and, more so, 
the ones that produce the most significant difference in 




Figure 15: Correlations between M 1/4 and E 1/12 . Top: The 
Signal for At = Mh = 120 GeV. Middle: WW background. 
Bottom: ti background. The horizontal scales are in GeV and 
do not have the same extent in all figures. 



16 




Figure 16: Top left: The distribution of E 1 / 12 for a signal with 
Mh = 120 GeV. The (yellow) histogram peaking at a smaller 
value of E is for the correct M — Mh- The other histogram 
is for M = 140 GeV. Top right: The distribution of M 1/4 
for "data" with M H 120 and M H 140 GeV. Bottom left: 
The distribution of E 1 2 for the same data, with M — Mh in 
each case. Bottom right: The same as the last entry, for Ai. 
In all cases the lower-peaking histogram his for Mh = 120 
GeV. The horizontal scales are in GeV. 



the shape of their distributions for a potential signal and 
the WW and tt backgrounds. 

Some of the singularity variables derived in the CM 
approximation of a motionless Higgs boson are unexpect- 
edly good. This is the case for Eo, associated with the 
CM condition of coplanarity, denned in Eq. (20) and il- 
lustrated in Figs. [5] and [6j Its two roots, M±, denned 
in the same equation, are of the simpler kind that does 
not involve a trial mass A4. Their correlations, shown 
in Figs. |7|8| in two different ways, are moderate, both for 
the signal and the background. The quantities M± and 
their product T%(M) = (M - M+){M - M_) are good. 

Still in the CM, we see in Fig. [9 that the naive vari- 
able Mi (or Ei), denned in Eqs. ( 22|24| ), is not good. 
The variable E2 (S3), also defined in Eq. ( |24| is good 
(not so good), as one can conclude from Figs. [9j [To] 
These limitations are lifted as we construct from these 
variables the quantities Di defined in Eq. ( [25] ) , which are 
histogrammed in Fig. [12] they are reasonably good at 
telling signal from backgrounds. Their correlation plots, 
shown in Fig. [l3j are quite disimilar for the signal and 
the WW irreducible background. 



The longitudinally boost-invariant analogs of D{ 




Figure 17: Study of the se nsit ivity to mass of A2, Eq. (36) 
(top row), and of A3, Eq. (37) (bottom row). The (yellow) 
histograms peaking at a smaller values are, in all graphs, for 
the correct M — Mh = 120 GeV. Left: Comparison with an 
incorrect choice M — 140 GeV. Right: Comparison with data 
for M H = 140 GeV, analized with the correct M = M H - The 
horizontal scale is in GeV. 



mass dimensionality in Eqs. (35 36|37 ). Only for Ai 5 2 we 
have analytical expressions, which are factorizable. The 
factors of Ai, and t he comple te variable are very good, 



14|l5|l6 So are the variables Ao 



17, 18 We see in Fig. 20 that the three 



as illustrated in Figs, 
as shown in Figs. 
Ai are quite correlated, but the correlation plots of signal 
and background are populated in a significantly different 
way. A comparison of Fig. [18] with the corresponding 
result for the CM variables Z^, Fig.[l2j shows, once more, 
that the A^ are demonstrably better. 

The confrontation of the results for the boost-invariant 
variables, A^, and their siblings, D^, obtained in the ap- 
proximation in which the boson is at rest is very gratify- 
ing. The signal peaks are significantly narrower and the 
correlations weaker for the As than for the Ds. 

In the sense of their correlations with the function Aip, 



shown in Fig. 21 , the singularity variables Ai are optimal 



1 to 3 are the singularity variables Ai of Eq. (29) and 
A2,3 of Eq. (31). We have redefined them to have unit 



tools to separate signal from backgrounds. 



XII. CONCLUSIONS 

Recall that, as discussed in the Introduction, in the 
case of the CM singularity variables one can construct up 
to five independent combinations of the relevant observ- 
ables. The best choice is the set {M + , M_, D i, £>2, D3}, 
whose ingredients are denned in Eqs. fl20|25| ). The mam 
appeal of these CM variables is that they are simple 
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Figure 18: Signal vs WW (left column) and it backgrounds 
(right column) for the variables Ai (top row), A 2 (mid- 
dle r ow) and A3 (bottom row), respectively defined in 
Eqs. fl35|36|37| ). In all graphs M = M H = 120 GeV. The 
horizontal scale is in GeV. 



and explicit functions of the relevant observables. Their 
main drawback is that they are not as good as the boost- 
invariant variables, to be revisited next. 

The only imperfection of the boost-invariant singular- 
ity variables is that for one of them, A3, we are unable to 
derive its explicit analytical expression. For A2 the an- 
alytical expression, Eq. (30), is so complex that we have 



opted to compute it event by event as a numerical resul- 
tant, as we are forced to do in the case of A3. For the 
computer, this is fast and simplest. 

The practical virtues of the variables Ai -in being able 
to pinpoint the actual value of M# and to tell apart signal 
from backgrounds- amply overcome their quoted single 
limitation. The theoretical toil required to go beyond the 
Higgs-at-rest approximation pays. 

Recall that in terms of the four boost-invariant observ- 
ables one can construct up to four useful combinations. 
The next-to-best choice is {M 1 / 4 , H 1 / 20 , A 2 , A 3 }, where 
M and £ are the factors building up Ai, see Eq. (29). 




Figure 19: High resolution histograms for A2 (upper panel) 
and A3, showing the (narrower, yellow) signal and the (wider, 
blue) WW background. The abscissae are in GeV units. 



The best choice is {Ai, A2, A3, Acp}, combining our kine- 
matical singularity variables with the good old "dynam- 
ical" (spin-dependent) angle, A(p, of the charged leptons 
in the transverse plane. 

The bell shapes of the signal histograms in Fig. 18 
are very satisfactory, even if obtained with theoretical 
expressions in the = approximation for the produced 
hadrons. We have checked, by generating and analizing 
events with = 0, that the improvement brought by 
theoretical variables that avoid the pt = approximation 
is unlikely to be very significant. 

We have only studied variables and their pair- wise cor- 
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Figure 20: Histograms of t he correlat ions between the singu- 
larity variables A t of Eqs. ( |35|36|37| ). Left: The M H = 120 
GeV signal. Right: WW background. Top: {Ai,A 2 }. Mid- 
dle: {Ai, A 3 }. Bottom: {A 2 , A 3 }. The A* axes are in GeV. 



Figure 21: Histograms of the correlati ons between the singu- 
larity variables Ai of Eqs. ( 35|36|37 ) and Acp. Left: The 
M H = 120 GeV signal. Right: WW background. Top: 
{Ai,A<p}. Middle: {A 2 ,A^}. Bottom: {A 3 ,A<p}. The A t 
axes are in GeV. 



relations. We have not attempted to quantify the abso- 
lute values of signals and backgrounds -as opposed to 
just the shape of their distributions. Thus, we are far 
from being able to show potential "significance" results 
in terms of a full multi-dimensional analysis of all vari- 
ables and their correlations. Yet our results for Ai 5 2,3 
in Fig. 18 are competitive in "goodness" with the Acp 



diagnosis recalled in Fig. [T] That was one of our goals. 

For an experimentalist eager to test the tantalizing 
hints that M H = 126 or 124 GeV DQ, it should not 
be too streneous to prepare the relevant one- or multi- 



dimensional singularity-variable templates for the rela- 
tively copious channel H — ^ WW —> leptons. 

Our main aim was the theoretical derivation of a com- 
plete set of phase-space singularity conditions and vari- 
ables for the process H -> W + W~ , W ± -» i^v. We 
have seen it is a rather laborious task. The origin of its 
difficulty is many- fold. First, because of the elusiveness 
of neutrinos, the kinematical constraints of Eqs. ([3| are 
incomplete. Second, several of these equations are non- 
linear. Finally and most severely, the 7-th equation, the 
one reflecting that the invariant mass of the four leptons 
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is Mh-, inextricably links the leptons resulting from the 
decay of one W to those from the other, very significantly 
complicating the ensuing algebra. 

In most processes relevant to a hadr on- collider search 
for new physics involving unobservable particles, the ini- 
tial step is a non-resonant production of a pair of novel 
particles. This means one cannot assume a fixed invari- 
ant mass for the pair and (approximately) boost each 
event to the pair's rest system. But the last difficultly 
mentioned in the previous paragraph is absent. That is 
why, even for Pt ^ -and a surfeit of unknown masses- 
the pertinent singularity variables are relatively simple, 
and analytical [TT] . 



Appendix: SO (1,1) invariance of the resultants 



Let d, Co be the CM functions defined in Eqs. fll8|19 ). 
Let C^Cq be the functions in Eq. (27), boosted by 
G SO(l,l), whose action on a longitudinal vector 
v = {^0,^3} is v ^ L^(v). Very explicitly, with £ the 
boost parameter in Eq. (128)), 



(38) 



v ^ (l 6 v) ■.= ^ + r>o + \(z - r>3, 
v 3 1 ^ (l 6 v) 3 ■.= 1(£ - r> + \(Z + r> 3 , 



and 



C' i (l,k,0 = CC i (L- 1 l,L^ 1 k), 
C' (l,k,O = C°C (L^l,L- 1 k), 



(39) 



where r^,ro are the minimal entire numbers required for 
C[, Cq to be polynomials in £. It is easy to check that 



i/2, r =g /2, 



(40) 



with gi and go the degrees in £ of C[ and Cq. 

Let R(l, k) be the resultant in £ of C[{1, k,£),C' (l,k,£): 

R(l, k) := Res(CK/, fc, 0, C^(Z, fc, 0,0- ( 41 ) 



We want to prove that 

= R(L v l,L v k) 



(42) 



for all G SO(l,l), that is, the resultant is invariant 
under longitudinal boosts. 

Indeed, using the aforementioned definitions, we have: 

R^Lrqlj L^k^j = 

Res{C[{L v l, L v k, 0, CS(L„i, L„k, 0,0 = 
ResCr'C^L" 1 ^/, L^L n k),C°C {L^L n l, L^L n k), = 
Res^C^L"/^, L-^k),e°C (L-^l, L^fc), = 

Res( v r <ci(i, k, z/ v ), v ro c^(i, k, s/ v ),z). 

(43) 

Taking into account that, if Pi(x),p2(^) are two arbi- 
trary polynomials of degrees g\ , #2 , respectively, 



Res(aipi(x/6), a 2 p2(x/b), x) 



af af 1 ft"* 1 * 2 Res(pi (a) , p 2 (x) , a) , 
it immediately follows, with use of Eqs. j40p4| , that: 
Resi^C'S, k,z/v), v ro ci(i, k, Z/t,), = 



(44) 



Res(Cf(l,fc,0,Ci(I,A!,0,0 = -R(i,A)- 



(45) 



The desired "Q.E.D ." is simply reached by putting to- 
gether Eqs. ( 43|45 ). It is also simple and gratifying, in 
the case of the A3 resultant that we were unable to derive 
explicitly, to check its boost invariance numerically. 
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